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Abstract 

We study chiral phase transition and confinement of matter fields in (2+l)-dimensional 
U(l) gauge theory of massless Dirac fermions and scalar bosons. The vanishing scalar boson 
mass, r = 0, defines a quantum critical point between the Higgs phase and the Coulomb 
phase. We consider only the critical point r = and the Coulomb phase with r > 0. The 
Dirac fermion acquires a dynamical mass when its flavor is less than certain critical value 
Nj, which depends quantitatively on the flavor Nb and the scalar boson mass r. When 
Nf < Nf, the matter fields carrying internal gauge charge are all confined if r ^ but 
are deconfined at the quantum critical point r = 0. The system has distinct low-energy 
elementary excitations at the critical point r = and in the Coulomb phase with r ^ 0. 
We calculate the specific heat and susceptibility of the system at r = and r ^ 0, which 
can help to detect the quantum critical point and to judge whether dynamical fermion mass 
generation takes place. 
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1. Introduction 

Quantum electrodynamics in (2+l)-dimensional space-time has wide applications in a 
number of planar strongly correlated systems, including fractional quantum Hall system, 
Heisenberg quantum antiferromagnet, and high-T c cuprate superconductor. Among its nu- 
merous variants, the massless spinor QED3 and the CP^ model are particularly interesting. 
They both exhibit fruitful physics in the context of quantum field theory and have proved 
to be powerful in understanding some fundamental physics of cuprate superconductors. 

It is widely accepted that the t-J model captures the essential physics of cuprate super- 
conductors [l|]. However, it is far from clear how to treat this model by satisfactory analytic 
tools. One frequently used approach to deal with the t-J model is based on the slave-particle 
technique [lj]. For example, the electron operator on z-lattice q can be decomposed as 

Ci = fA, (1) 
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where fa represents the operator of spin-carrying spinon and b\ represents the charge-carrying 
holon. Such decomposition has a redundancy since the electron operator is invariant under 
the transformations 

fi^e i0 f h h^e ie b h (2) 

which can be considered as an emergent gauge symmetry. This emergent gauge structure has 
very important effects on the whole physical picture and has been investigated extensively 
for twenty years [H, Q, 0, 0, H, 0]- However, such decomposition is not unique and one can 
alternatively use the so-called slave-fermion approach, in which the spin quantum number 
is carried by scalar field [tJ and the gauge charge is carried by fermionic matter field [sj]. 

Recently, a new state of matter, named as algebraic charge liquid, was proposed to 
describe the low-energy properties of cuprate superconductors 0, 0]. The effective field 
theory of this state consists of iVy-flavor massless Dirac fermions, iVVflavor scalar bosons, 
and a U(l) gauge field. This effective theory can be loosely considered as a combination 
of massless spinor QED 3 and CP^ 1 model 0|. It has been used to understand certain 
quantum phase transitions and various exotic properties near the quantum critical points 
ji], 0]. Recently, this model was studied by Kaul and Sachdev [oj with emphasis put on the 
anomalous dimensions and some observable thermodynamic quantities including specific 
heat and susceptibility near the quantum critical point between Higgs and Coulomb phases. 

When interacting with U(l) gauge field, there is an interesting possibility that the mass- 
less Dirac fermion mig ht undergo a vacuum condensation, which then generates a finite mass 
for Dirac fermion [loL ElL 12]. This mass breaks the continuous chiral symmetry possessed 
by Dirac fermions and leads to new phenomena. Another characteristic property of QED3 
is that the strong gauge field can confine the matter fields that carry nonzero gauge charges 



13l . 1 14T | . In the absence of scalar field, it is found that the fermion mass can be generated 
only for flavor Nj < NS « 3.3 [10], and confinement takes place once the mass is generated 



141 ] . If the gauge field also couples to the scalar field z, then the phenomena related to chiral 
phase transition and confinement will be modified by the dynamics of scalar bosons. 

In this paper, we study chiral phase transition and confinement of matter fields in the 
vicinity of the quantum critical point. The critical fermion flavor Nj depends on two pa- 
rameters: the scalar boson flavor Nj, and the scalar boson mass r. When r > 0, the system 
stays in the Coulomb phase, while r = defines the quantum critical point that separates 
the Coulomb phase from the Higgs phase. When the chiral symmetry is dynamically broken, 
the massless fermions become massive and are confined if r 7^ 0; while the Dirac fermions 
and z bosons are deconfined at the quantum critical point r = 0. Since the confined and 
deconfined matter fields behave quite differently, the Coulomb phase and the critical point 
have completely different low-energy elementary excitations. Thus it is possible to specify 
the quantum critical point and to verify whether the chiral symmetry is dynamically broken 
by measuring observable physical quantities such as specific heat and susceptibility. We 
calculate these quantities within the effective gauge theory and discuss their behaviors at 
low temperature. 

The paper is organized as follows. We study the dynamical mass generation for Dirac 
fermions in Section 2. The confinement /deconfinement properties of matter fields is dis- 
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cussed in Section 3. We then calculate the specific heat and susceptibility at quantum 
critical point r = in Section 4 and in the Coulomb phase with r ^ in Section 5. The last 
section contains summary and discussion. 



2. Dynamical fermion mass generation for r > 

The action for the interaction between massless Dirac fermion and U(l) gauge field is 



S f = \f d 3 xF^ + Y,f d 3 x^ % - 



ia^ipi, (3) 

where the fermion has Nf flavor. The spinor field ip describes the superconducting quasi- 
particle that carries physical electric charge. Here, for generality we keep the kinetic term 
of gauge field. Following Kaul and Sachdev we consider only non- compact gauge field 
and omit the instanton excitations, so the above field theoretic model might not precisely 
match the low-lying theory of cuprate superconductors (there are theoretical and simulation 
evidence for the irrelevance of instantons in the presence of massless Dirac fermions with 



sufficiently large flavor Nf [15l . I la , Il7fl). In order to define chiral symmetry, we adopt a 
four-component representation of spinor field ip and three 4x4 matrices 7^ (fi = 0,1,2) 
satisfying the Clifford algebra, {7^,7^} = 25 M „. Since the fermions are massless, the action 
Sf respects the following continuous chiral symmetry 

i) -> e^^ip, (4) 

with 73 and 75 being two 4x4 matrices that anticommute with 7^. The action contains no 
small coupling constant and generally can only be treated using the 1/Nf expansion. 
The free propagator for Dirac fermion is 

G^\p)=ii^. (5) 
The interaction with gauge field modifies it into the complete fermion Green function 

G-»=z7^ + £(p), (6) 

where E(p) is the fermion self-energy function. Generically, the fermion self-energy can be 
written as 

E(p) = i^(A(p 2 )-l) + B(p 2 ), (7) 

with A(p 2 ) being the wave functional renormalization and B(p 2 ) the fermion mass function. 
To calculate the fermion self-energy, we can use either conventional perturbation expansion 
method or non-perturbative approach. If the usual perturbation theory is used, then the self- 
energy function is proportional to ~ 777 -pln(p) to the leading order of 1/Nf expansion 0, 0], 
which represents only a correction to the wave function (field operator). When combined 
with the free propagator «7^p M , this self-energy modifies the fermion propagator to a new 
form 

G{p) - 
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which reflects the typical behavior of algebraic charge/spin liquid [y, This result implies 
that, to the leading order, the gauge fluctuation is only a marginal perturbation and no mass 
term B(p 2 ) is generated. In fact, such mass term can never be generated to any finite order 
of perturbative expansion. This can be understood by the symmetry arguments. Both the 
free action and the term i r j^p^A{p 2 ) are chiral symmetric, but the mass term B(p 2 ) is 
not chiral symmetric. Therefore, it is the chiral symmetry of the action that prevents the 
occurrence of mass term B(p 2 ). 




Figure 1: The solid line represents the propagator of Dirac fermions, dashed line the scalar boson z, and 
dotted line the field A. (a) and (b) correspond to the polarization function from z bosons, (c) corresponds 
to the polarization function from Dirac fermions. (d) corresponds to the polarization of field A. 



However, the perturbation theory can not tell us everything. Although the action is 
chiral symmetric, it is possible that its ground state does not respect the chiral symmetry 
and the Dirac fermion acquires a finite dynamical mass. Indeed, the generation of nonzero 
B(p 2 ) signals the happening of a phase transition and the reconstruction of the ground 
state, which are surely non-perturbative in nature and can not be analyzed by ordinary 
perturbation theory. The standard non-perturbative approach to the problem of dynamical 
fermion mass generation is to solve the Dyson-Schwinger (DS) equation for the fermion 
mass function, which was originally proposed by Nambu and Jona-Lasinio in the context of 



particle physics [18|. The DS equation can be written in the form 



^G(k)r v (p,k)D„ u ( P -k), 

where T v {jp, k) is the vertex function and D^ v {p — k) the photon propagator. After a series 
of manipulations, the above equation is decomposed into a couple of equations for A{p 2 ) 
and B(p 2 ) respectively 

73/ 



MP 2 ) = ^ + ^j^f^b-kl ll G{k)V v D lxv {p-k)l 
B(P 2 ) = -~j ^Tr^Gik^D^p-k)}. 



If the DS equation for B(p 2 ) has only vanishing solution, the fermions remain massless 

and are stable against gauge fluctuations. Whenever the equation develops a nontrivial 
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solution B(p 2 ), the massless fermion acquires a finite mass. The crucial difference between 
the perturbative and non-perturbative approaches is that the DS equation is nonlinear which 
makes it possible for chiral phase transition to happen at the bifurcation point. 

The logic presented above is analogous to the gap generation for fermion excitations in 
the BCS theory of conventional superconductors. In a normal metal, when we compute the 
mean value (c^c^ki) using ordinary perturbative expansion, only vanishing result can be 
obtained. It could be understood that the local gauge symmetry of the system prevents the 
appearance of nonzero (c^C-ki) which is not gauge invariant. However, a finite gap can 
be obtained by solving the nonlinear gap equation that is constructed by non-perturbative 
approach. Apparently, the local gauge symmetry does not prevent the gap generation; on 
the contrary, this symmetry is spontaneously broken by the finite gap since (c^c^ki) 7^ 0. 
The gap generation signals a phase transition from normal metal to superconductor. In fact, 
the idea of dynamical fermion mass generation was motivated by the Cooper pairing picture 



in BCS theory [18j. The order parameter for chiral phase transition is just the fermion mass 
term (ipip) which is caused by fermion vacuum condensation. 

For the spinor QED model ([3D, the chiral phase transition has been studied for many 
years. It is known that there exists a critical fermion flavor Nf which separates the chiral 
symmetric phase (Nf > Nf) from the chiral symmetry broken (Nf < Nf) pha se < 10l . I 111 Il2 



The critical flavor is about Nj ~ 3.3 to the leading order of 1/Nf expansion [1( 

We now add the action of complex scalar field. At half-filling, the cuprate supercon- 
ductor is a Mott insulator with long-range Neel order. The Neel order parameter n can be 
represented by the scalar boson operator z a by 

n = z*a a pZfs. (9) 

Then the Neel state can be described by an effective field theory called CP^ 6-1 model, where 
Nb is the flavor of scalar field z^. This model admits an emergent local gauge symmetry, 



as demonstrated by Polyakov [19|, and has been successfully used to describe the planar 
Heisenberg quantum antiferromagnetism @, @, 0] • Its action is 0, @, 0] 

Sb = ~ g J (Px^-ia^ 2 -i\(\z a \ 2 -1)] . (10) 

The Lagrangian multiplier A is introduced to impose the constraint \z a \ 2 = 1. It takes on a 
uniform saddle point value that extremizes the action A = ir, 

d3 » 1 1 no 



(27r) 3 p 2 + r g 



For g > g c , the system stays in the Coulomb phase with (z a ) = 0; for g < g c , the system is in 
the Neel phase where the SU(Nfe) rotational symmetry is broken by (z a ) 7^ 0. In this paper, 
we only consider the phase with (z a ) = 0. At the Nf, — > 00 limit, the quantum critical point 
locates at r = 0, where g = g c satisfying 

d3 » 1 1 (12) 



(2vr) 3 jo 2 g c 
5 



From these formula, we have 



l_ 1 f d 3 p ( 1 1 \ y c? 3 p r _ y/r 

# c g~ J (2tt) 3 VP 2 p 2 + ry 7 (2vr) V(p 2 + r) ~ 4vr ' 1 j 

from which the mass parameter 

A 2 (9 - 9c? 



7r 2 g 2 



(14) 



For finite Nf,, the scalar field z receives self-energy correction from gauge fluctuations, so the 
critical point is no longer at r = 0. Even in this case, one can define a new quantity r g so 
that the system goes critical as it vanishes ||. Here, we simply use the single parameter r 
to denote the boson mass. 

The additional coupling between gauge field and scalar field will affect the critical fla- 
vor NJ. The DS equation of fermion self-energy function will be studied within the 1/Nf 
expansion. In the Landau gauge, the gauge boson propagator has the expression 

D - (g) = g 2 [l + n a (g 2 )](^-^)' (15) 
where the polarization function is given by 

n a (g 2 ) = n{(g 2 ) + n^ 2 ), (16) 



with the fermion contribution Il{(g 2 ) = Nf/8\q\ (Fig. 1(c) ). The scalar sector will be studied 



using the l/iV& expansion. In particular, the scalar boson contribution to one- loop vacuum 



polarization (Fig. 1(a) and Fig. 1(b)) is 



Inserting the propagator D^ v (q) into the mass equation, then the critical fermion flavor Nf 
becomes a function of r and Nf,. We now determine this dependence by solving the mass 
equation. 

To the lowest-order, the wave function renormalization can be simply taken to be A(p 2 ) = 
1 since the next-leading order is suppressed by a factor of 1/Nf. In the Landau gauge, the 
vertex function reduces to Y^{p,k) = 7^, as required by the Ward identity. Under these 
approximations, the mass equation is relatively simple. To go beyond this approximation, 
one should solve the coupled equations for A{p 2 ) and B{p 2 ). Further, once A{p 2 ) ^ 1, 
the vertex function can no longer be taken to be 7^. It must be chosen properly to be 
consistent with the Ward identity. As it turns out, the DS equations become much more 
complicated after including these higher order corrections. It is also difficult to chose a 



proper vertex function. Fortunately, there is a very useful technique [2Cj, [2JJ which simplifies 
the calculation. To use this technique, we first write the following propagator for gauge 
boson 

O 



where £(<? 2 ) is a gauge fixing parameter that depends on momentum. Taking advantage 
of the gauge symmetry, one can choose a so-called non-local gauge £(<? 2 ) so that the wave 
function renormalization is strictly unity, A(p 2 ) = 1. Then the vertex function becomes 
^ flip j k) = j p. The gauge parameter is determined by the following expression 

Now we only need to solve a single DS equation for mass function 

B(n*) f d " k ^ 3-^(( P -A;) 2 ) 
KP ' J {27r) 3 k 2 + B 2 {k 2 ){p-k) 2 {l + U a {{p-k) 2 ))' y ' 

This integral equation is much more easily solved numerically than the coupled equations 
of A{p 2 ) and B{p 2 ). It reduces to the mass equation in Landau gauge when the gauge 
parameter £ = 1. There must be a critical fermion flavor Nf below which a finite fermion 
mass is generated. 

We solved the DS equation j22| in both the Landau gauge and the non-local gauge, with 



results presented in Fig. 2(a) and Fig. 2(b) respectively. Note that the parameter r now 
becomes dimensionless since it is re-scaled by some intrinsic energy scale (the parameter a 
defined in the paper of Appelquist et al. [lOj). The results show that the boson mass r 
increases the critical fermion flavor N* and that the boson flavor Nb reduces Nf for small 
mass r. When r is of the order unity, the effect of scalar boson z a on Nf drops out and the 
critical flavor Nf is independent of r and Nf,. This is easy to understand because for large 
r the scalar boson contribution to polarization LT^ — > N^/r can be neglected compared with 
the fermion contribution Il{(g 2 ) = Nf/8\q\ at low momenta. If we believe that the results 
obtained in the non-local gauge is more reliable than that in Landau gauge, then the critical 
flavor Nf > 2 for iVj, < 4. When the physical fermion flavor Nf = 2, a finite fermion mass 
is always generated for Nf, < 4. 

Once the fermion mass is generated, the system is no longer critical even at the quantum 
critical point r = 0. The nonzero value of (ipip) corresponds to the formation of charge 
density wave since the fermion field ip carries only electric charge. As the result of dynamical 
breaking of continuous chiral symmetry, there appear massless Goldstone bosons which 
are fermion-antifermion pairs. The finite fermion mass has substantial effects on various 
observable physical quantities, to which we turn in the next three sections. 

When the scalar bosons condense, (z a ) ^ 0, the gauge boson becomes massive due to 
Anderson-Higgs mechanism. The effect of gauge boson mass on has been discussed 



previously in Ref . [22] , where it was found that this mass rapidly suppresses the dynamical 
mass generation for Dirac fermions. We will not repeat such considerations in the present 
paper. 

When applying the effective field theory to realistic cuprate superconductors, the Dirac 

fermions actually have two velocities, vp and v a, with vpjv^ ~ 20 [lj. Although the velocity 

anisotropy is important in certain physical quantities, it was shown [23| to be irrelevant near 

the critical point of chiral phase transition. Indeed, in jnost^ubhcations on spinor QED 3 

the velocities are simply set to unity, i.e., vp = va 
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in most publica 




(a) (b) 

Figure 2: (a) Dependence of critical flavor N'j on scalar boson mass r for different values of boson flavor 
N\, in the Landau gauge, (b) Dependence of critical flavor on scalar boson mass r for different values of 
boson flavor in the non-local gauge. 



3. Classical potential of gauge field: confinement and deconfinement 

The matter fields interact with each other through the mediation of U(l) gauge field. 
The classical potential between two particles with opposite gauge charges can be written in 
the coordinate space as [13] 

v ^=-lw? e "w^mx (21) 



It has the following asymptotic form [13| 

V(x)~- L— Wlxl). (22) 

v ; 27ri + n(o) u u K J 

Obviously, the potential is determined by the infrared behavior of vacuum polarization 
function. In the absence of vacuum polarization II(g 2 ), the potential has a logarithmical 
form 

V(x) ~ ln(|x|). (23) 

This potential increases at large distances and hence is a confining potential, though weaker 
than a linear potential. When the polarization II(g 2 ) is included, the gauge potential will 
be changed. 

In cases without scalar field, N b = 0, the form of gauge potential V(x) depends only on 
whether the fermion mass is zero or finite 14| . For massless fermion, the vacuum polarization 
is II(g 2 ) = gj4. Since 11(0) — > oo, the confining potential is destroyed. After the fermion 
mass is generated, the polarization function Il{(g 2 ) becomes 

Nf ( q 2 - 4m 2 

' ' arcsm ■ 



n£( « 2) = 4^( 2m + 



4vrg 2 V q ~ i/o 2 + 4m 2 

8 V 



In the limit q — > 0, it has a finite value 



n{(0) = ^. (24) 
Thus in the symmetry breaking phase, the gauge potential is confining and there can not 



be any asymptotic fermion states 14 



Note the gauge field considered here is non- compact and there are no topological in- 
stanton configurations. The confinement is induced completely by the gauge force between 
particles with opposite charges. It is a special feature of gauge field in two spatial dimensions. 

However, the gauge potential V(x) is modified by the additional scalar boson. Depending 
on the fermion mass m and the scalar boson mass r, there are four possibilities: 

1. if m = r = 0, II(g 2 ) = 2N 1 f e ^ b diverges at q — > 0; 

2. if m ^ and r = 0, 11(0) = 5 + ifu diverges at q -> 0; 

3. if m — and r^0, 11(0) = gj4 + ^ diverges at g — ► 0; 

m ' r 



4. if m ^ and r ^ 0, n(0) = ^ + ^ is finite. 



Compared with the case without scalar bosons, there appears a new possibility: at ex- 
actly the critical point where r = 0, the gauge potential is always deconfined, irrespective 
of whether the chiral symmetry is broken or not. This implies that, in the chiral symmetry 
broken phase, the matter fields carrying internal gauge charges are deconfined at the quan- 
tum critical point r = 0, but are always confined once the system enters the Coulomb phase 
with r ^ 0. The deconfined matter fields will make finite contribution to observable physical 
quantities, such as specific heat and susceptibility, but the confined matter fields make no 
contributions to them. In the confining phase, only excitations with zero gauge charge can 
be mobile and make finite contributions to observable quantities. Obviously, the quantum 
critical point and the Coulomb phase have rather different low-energy elementary excita- 
tions. It is possible to distinguish the critical point from the Coulomb phase by measuring 
the observable quantities. These quantities can also be used to judge whether dynamical 
chiral symmetry breaking happens or not. 



4. Deconfinement and observable quantities at the quantum critical point r = 

We first discuss the quantum critical point, r = 0, where the Dirac fermions and z 
bosons are deconfined. In this point, the low-energy elementary excitations are: massive 
Dirac fermion if>, massless scalar boson z, and composite massless Goldstone boson. The 
gauge field is strongly interacting with Dirac fermions and scalar bosons, so it makes finite 
contributions to the specific heat and susceptibility. 

We first calculate the free energy T of the system, which is given by 

F = N b f b + N f ff + / G , (25) 

where f ob , f°f, and f G are the free energy of free scalar bosons, free fermions, and Goldstone 

bosons, respectively. To calculate the susceptibility, we follow the strategy of Ref. [i| 

9 



and introduce two magnetic fields, H b and Hf. For z bosons, the field shifts frequency 
as u n — > u n — 6H b /2, where 9 = ±1. For fermions, the coupling is the same, with the 
replacement u n — > u n —9Hf. Then the susceptibility x an d specific heat Cy can be computed 
by 

Xb = dH~r (26) 
* = as* (27) 

<9 2 .F 

Cv = -T^r 2 . (28) 

If the system is completely critical, the specific heat has the typical behavior ~ T 2 . Once 
the Dirac fermion gets massive, the system is surely no longer critical and the ~ T 2 behavior 
must be changed. In order to make comparison, we still write the specific heat at low-T in 
the form 

C v = A Cv T\ (29) 



A Cv = N b A° c b v + N f A% +AZ, (30) 



where the coefficient 

'%+N f A°l+A G Cv 

depends on temperature T. The susceptibility is linear in temperature T when the system 
is critical, but also deviates from this behavior when the fermion becomes massive. It has 
the form 

= A%T, (31) 

for scalar boson z and 

Xf = A%T, (32) 

for fermion ip. The susceptibility for Goldstone bosons certainly vanishes since they do not 
couple to any magnetic field. 

4-1- z scalar bosons 

When the system is not at the critical point, r ^ 0, all matter fields carrying internal 
gauge charges are confined. Therefore, their contributions to specific heat Cy and suscep- 
tibility x° b strictly vanish at low-T. We only need to consider the point of r = 0, where the 
free energy is given by 



f° b = | E / Sm^ 2 + K + H b /2f + ml) + ln^ 2 + (u n - H b /2) 2 + m 2 b )] - 
The parameter m b depends on T and field H b as 0] 
m b = Tin 



mi 



(33) 



2 cos{H b /2T) + 1 + v /{2cos{H b /2T) + l) 2 - 4 
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The contributions of z bosons to specific heat and susceptibility have been computed by 
Kaul and Sachdev with the results 

A°S V = 6^ « 1.83661, (34) 

j& = A n ^[±I~ 0.0856271. (35) 

4-2. Massive fermions 

At low but finite temperature T, the fermion mass should be a function m(T, k) of T 
and momenta k. It is usually difficult to get an analytical expression of m(T, k). As an 
approximation, we replace this function by a constant m, which might be identified by the 
value m(0, 0) at the limits T = 0, k = 0. When r = 0, the free energy of massive fermions is 

ff = _ T /^[Wi + e -(VF+^+iW^ (36) 

J 4:TC 2 

from which the specific heat of massive fermions is found to be 



A of - - 



^) 3 — ^ + 3(-) 2 ln[l + eT] + 6(-)Li 2 [-eT]-6Li 3 [-eT] . (37) 
1 1 + eT 1 1 



Here, Li n [2] is polylogarithm function. 
At the limit m — > 0, it reduces to 



_5 Li3 [_ 1]=6 ?g). (38 ) 

71 47T 



The dependence of A°q on ^ is depicted in Fig. 3(a) The susceptibility of massive 
fermions is 

i m — 

A% = V^+Ml + e*]). (39) 
As m — > 0, it approaches — . For finite m, the dependence of ^4.^ on ^ is shown in Fig. 



3(b) It is obvious that the mass suppresses the fermion contributions to both specific heat 



and susceptibility. Comparing with the contribution to specific heat from the z bosons, the 
contribution of massive fermions can be simply neglected. 

4-3. Goldstone bosons 

Since the composite Goldstone bosons are neutral with respect to gauge charges, they are 
not confined by the gauge potential and can make finite contributions to certain observable 
quantities. At r ^ 0, the Goldstone bosons are the only degrees of freedom at very low- 
energy regime. Their contributions to specific heat are nonzero and have the same expression 
at both r = and r / 0. The free energy for Goldstone bosons is 

/ G = r £/0H* 2 +^,)] = -^. («) 
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(a) Specific heat of massive fermions at very low T . The top dashed line corresponds to J^q = 
Susceptibility of massive fermions at very low T. The top dashed line corresponds to A^f = 



The specific heat is A% y T 2 with 



A, 
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C(3) 



4.591519. 



(41) 



T dT 2 7T 

Their contribution to susceptibility is zero because they do not couple to any magnetic field. 

4-4- Gauge contributions to physical quantities at critical point r = 

At the deconfined critical point, r = 0, the gauge field might give rise to important 
contributions to specific heat and susceptibility. The effective theory contains both temporal 
and spatial components of the gauge field, so there is no singular corrections to the free 
energy [1, [24| . 

To calculate the gauge field corrections to free energy, we should first obtain the polar- 
ization functions. The Feynman diagrams for the fermion and scalar boson contributions 
to polarization functions of gauge field are presented in Fig. [TJ At finite temperature, the 
one-loop vacuum polarization tensor are defined as 



n L(Po,p) 



nL(po,p) 



n 

El V 



d 2 k TrfrJ-yM- 



(2tt) 2 

d 2 k 

(2^)2 



k 2 (p + k) 2 ' 

(p + 2fc) /1 (p + 2fc) t 
k 2 (p + k) 2 



s, 



~k? 



(42) 



Here, the Matsubara frequency is p = 2rmr/ j3 for gauge boson, k = 2nn / ' (3 for scalar boson, 
and the frequency is k$ = (2n + 1)tt/(3 for Dirac fermion. Using the transverse condition, 
p^H^u = 0, the polarization tensor can be decomposed as 



lWpo,p) = [U{(p ,p) + U b A (p ,p)}A^+[U f B (p ,p) + U b B (p ,p)}B llu , (43) 
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where 



A - f X PfPo\ P 2 f X P°P» 



Bp, = 5 lxi [8 ij - V ^-)8 jv , (44) 



which satisfy the relationship 



V + B , v = <W " Vj j-- (45) 



Now the gauge boson propagator can be written as 

A«,(po,p) = f — + f — • ( 46 ) 

+ n^( Po , P ) + u b A ( Po , P ) p> + n / (p0j p) + n ^ (po; p) 

The gauge contribution to free energy is given by: 

/ 1A (A H f ) = T - J ^ [(P 2 + Hi + IW + 4 + 11$,)] . (47) 



£ m 



The vacuum functions 11^, IT^, IT^, and IT^ should be calculated explicitly They are related 
to the temporal and spatial components of vacuum polarization tensor li.^ and H b by the 
identities 

n$» = gnff\ (48) 
ng » = nf»-^n£" (49) 

In Ref. 0, the temporal and spatial components of polarization functions were obtained 
approximately. Indeed, these functions can be computed exactly even in the presence of 
finite magnetic field Hf^ and fermion mass m. Using the methods presented in Ref. 251 ]. 



we obtained the following expressions for the fermion contribution to polarization functions 

l4( em ,k) = *l f dx\n[4D m (x,e m ,k,H f )]-^L f dx{1 _ 2x)e J^f €m - 

ii(3 J 4vr J D m (x,e m ,k,H f ) 

N f f 1 , m 2 + x(l-x)ei , . 

/ dx - - , -+-g-swh(pE f ), 

27tJ q E f D m {x,e m ,k,H f ) KH th 



H^^m; k) 



_Nf_ j l dx ()_-2x)e m sm(x0e m - 09H f ) 
Jo 



D m {x, 6 m , k, H 



f) 
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N f [\ x(l-x)(k 2 + 2e 2 m ) . . 
Att J E f D m (x,e m ,k,H f ) KH i! 



The corresponding expressions for scalar bosons are 

ll 00 (e m ,k) = — — / cfeln [AJ m (x, e m ,k, H b )\ - — / da;(l-x)e r 
Jo 47r 7o 



Jmk^i €mi k, -Hft) 



4tt J E b J m (x, e m , k, f/ 6 ) 



n-ii(e m , k) — n 22 (e m , k) 

N b f 1 , xe m sm(x/3e m -§8H b ) N b f 1 x(l - x)(k 2 + e 2 J + \k\ - x 2 e 2 m . 

dx — : — — h— / dx — — smh(/lE; 6 ). 



47T J J m ( 

2^; £tnj k, Hh) 8ll Jq E b J m {x i 6 m , 

Here, we have defined a list of functions: 

D m (x, e m , k, H f ) = cosh 2 {^(3E f ) - sm 2 (^xf3e m - ^0H f ), 

E f (x,e m ,k,m) = ^ 'm 2 + x(l - x)(k 2 + e£J, 

1 1/3 
J m (x, e m , k, = cosh 2 (-/3£ 6 ) - sin 2 (-x/3e m - -9H b ), 

E b (x, e m , k) = y/x(l - x){k 2 + e|J. 

Note the scalar mass is taken to be r = since the gauge field is deconfined only at this 
critical point. 

The field A has dynamics only by integrating out the boson fields z, which leads to the 
following free energy 

f lx W,H b ) = ^E/(0i ln tnA(6 m ,k)], (so) 



where the vacuum polarization, given by Fig. 1(d) , at finite temperature is 

16vr Jo E b J m (x, e m , k, H b ) 

Now we are ready to calculate the gauge field corrections to free energy f lA and f ix 

by inserting the polarization functions H^, 11^, 11^, 11^ into (14"T|) and U\ into (I5D1 . Then 

the specific heat and susceptibility can be obtained by taking derivatives of the free energy 

with respect to T and Hf b , respectively. After numerical computations, the results for 

specific heat and susceptibility are presented in Fig. H] and Fig. [51 Apparently, the gauge 

contributions to specific heat and susceptibility are significantly suppressed by the mass 

of Dirac fermions. Here, in order to show the explicit results, we use the absolute value 

of \A$ V \ = -A l ^ v (\A£\ = -A 1 ^) instead of A£ v (A£). Although r = at the critical 

point, for completeness we also present the results for some magnitudes of finite r in Fig. 

[51 If Nf > N*, the fermions remain massless and the matter fields carrying internal gauge 
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charges are not confined at finite r. In such state, the mass r suppresses the specific heat 
and susceptibility of z bosons. 

In summary, at the critical point r = 0, both fermionc and gauge field contributions 
to specific heat and susceptibility are suppressed by the dynamical fermion mass. At low 
temperature, the total specific heat receives contributions mainly from three parts: free 
scalar boson z, the composite Goldstone boson, and the fluctuation of field A. Since the 
z and Goldstone bosons are both massless, these three contributions to specific heat are 
all proportional to T 2 , with coefficients Ac v being temperature independent. However, the 
scalar bosons carry nonzero internal gauge charge, while the Goldstone bosons are neutral. 
Therefore, the former lead to a oc T susceptibility (to H b ), while the latter does not couple 
to any magnetic field {Hf, Hf,, or the physical H) and hence make no contribution to 
susceptibility. 




4>) 




(a) 



(b) 



Figure 4: (a) The gauge field contribution to specific heat for different fermion mass, (b) The gauge field 
contribution to susceptibility for different fermion mass. 




10* 10" e T 10 -1 10" z 10"* 10* T TO - ' ID -2 



(a) (b) 



Figure 5: (a) The contribution of field A to specific heat for different boson mass, (b) The contribution of 
field A to susceptibility for different boson mass. 
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5. Confinement and observable quantities at r 7^ 

We now consider the Coulomb phase with r / 0. Once the logarithmic gauge potential 
confines all particles that carry internal gauge charges, the excitations appearing at low 
energy regime must be neutral with respect to internal gauge charge. Now the massive 
fermions can combine either with anti-fermions or with massive scalar bosons carrying op- 
posite gauge charges. In the former case, the corresponding composite particles are just the 
massless Goldstone bosons. In the later case, the separated spin and electric charge degrees 
of freedom are re-combined by the confining internal gauge potential. This implies that 
there are also electron-like quasiparticles that carry both 1/2-spin and — e electric charge. 
It is a nontrivial task to estimate their mass since this requires a clear understanding about 
the mechanism of spin-charge recombination. The issue of s pin -charge recombination has 
been investigated for many years by several authors (a, IH, E5, l27||. Unfortunately, this issue 
is not at all clear. In the absence of a reliable theory for spin-charge recombination, here 
we assume a constant mass m e . Further, the z bosons can combine with the z* bosons to 
form composite spin excitations 0, @, B). Both the electron-like fermionic quasiparticles 
and the composite spin excitations are gapped, so their contributions to specific heat and 
susceptibility (with respect to real magnetic field H, rather than £//,&) are suppressed by 
their mass gaps m e and A spin . At very low temperatures, only the composite Goldstone 
bosons contribute to specific heat (see equation (41)). However, although having the same 
T-dependence at the point r = 0, the magnitude of specific heat in the Coulomb phase is 
much smaller than that at r = since the contributions from massless z bosons and the 
field A all disappear at r 7^ due to confinement. 

6. Summary and discussion 

In summary, we have studied the possible mass generation of initially massless Dirac 
fermions in a continuum U(l) gauge field theory in (2+1) dimensions. It was found that a 
dynamical mass is generated when the fermion flavor Nf is less than which has quantita- 
tively different values at the quantum critical point r = and at r 7^ 0. The gauge potential 
between charged particles leads to logarithmic confinement when r ^ 0, but at the critical 
point r = the matter fields carrying nonzero internal gauge charge are deconfined. The 
low-energy properties of the system at r = are quite different from those at r 7^ 0. We dis- 
cussed some observable physical quantities, including specific heat and susceptibility, of the 
system at low temperature and argued that they can be used to judge whether the system 
stays at the quantum critical point r = and whether the fermions acquire a dynamical 
mass. 

Besides specific heat and susceptibility, the happening of dynamical fermion gap gen- 
eration can also be readily seen from the transport properties of the system. At very low 
temperature, the thermal conductivity is mainly determined by the low-lying fermionic ex- 
citations. If the Dirac fermions are massless, then they exhibit universal thermal conductiv- 
ity, ~ = I V ^fva ' wmcn i s independent of impurity scattering rate r jig}. Once the Dirac 
fermions acquire a dynamical mass m, then the thermal conductivity at zero temperature 
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becomes & = 4 — — - „■> ,° 9 , which is no longer universal and is reduced by the fermion mass 



29[. If the system is in the clean limit so that T <C m, then the thermal conductivity 
vanishes at very low temperature. 

When chiral symmetry is dynamically broken (no matter r = or r / 0), there are 
no gapless fermionic excitations at low temperature. Therefore, in the limit ro ■€ m the 
thermal conductivity effectively vanishes at zero temperature and does not have a universal 
linear term that would exist if the system contains gapless Dirac fermions. The ground 
state of the system corresponds to an interesting state of matter where the specific heat ex- 
hibits behavior of gapless (bosonic) excitations but the thermal conductivity has no residual 
linear term. In two recent works, thermodynamic measurements carried out in k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 found evidence for the presence of gapless excitations [30], while thermal 
transport measurements carried out in the same material found good evidence supporting 



the absence of gapless fermionic excitations 31]. Although it is not known whether the 



ground state of k-(BEDT-TTF) 2 Cu2(CN) 3 can be described by a continuum field theory 
analogous to the combined action Sf + Sb, such seemingly contradicting experimental facts 
are consistently realized by the chiral symmetry breaking phase of action Sf + Sb and can 
be naturally understood by the mechanism discussed in the context. 

The analysis presented above can also be applied to the Abelian Higgs model, which may 
be considered as the relativistic version of Ginzburg-Landau model, with Higgs field being 
the superconducting order parameter. Although the CP Nb ^ 1 model is distinct in physical 
contents from the Abelian Higgs model, its low-energy property bears much resemblance to 
that of the latter. 



Acknowledgments 

This work is supported by National Science Foundation of China No. 10674122. 



References 

[1] P.A. Lee, N. Nagaosa, X.-G. Wen, Rev. Mod. Phys. 78 (2006) 17. 

[2] I. Affleck, J.B. Marston, Phys. Rev. B 37 (1988) 3774. 

[3] L. B. Ioffe, A.I. Larkin, Phys. Rev. B 39 (1989) 8988. 

[4] D.H. Kim, P.A. Lee, X.-G. Wen, Phys. Rev. Lett. 79 (1997) 2109. 

[5] D.H. Kim, P.A. Lee, Ann. Phys. (N.Y.) 272 (1999) 130. 

[6] W. Rantner, X.-G. Wen, Phys. Rev. Lett. 86 (2001) 3871. 

[7] T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, M.P.A. Fisher, Science 303 (2004) 1490; T. Senthil, 

L. Balents, S. Sachdev, A. Vishwanath, M.P.A. Fisher, Phys. Rev. B 70 (2004) 144407. 
[8] R.K. Kaul, A. Kolezhuk, M. Levin, S. Sachdev, S. Senthil, Phys. Rev. B 75 (2007) 235122; R. K. Kaul, 

Y.B. Kim, S. Sachdev, S. Senthil, Nature Physics 4 (2008) 28. 
[9] R.K. Kaul, S. Sachdev, Phys. Rev. B 77 (2008) 155105. 
[10] T. Appelquist, D. Nash, L.C.R. Wijewardhana, Phys. Rev. Lett. 60 (1988) 2575. 
[11] D. Nash, Phys. Rev. Lett. 62 (1989) 3024. 

[12] E. Dagotto, J.B. Kogut, A. Kocic, Phys. Rev. Lett. 62 (1989) 1083. 
[13] C.J. Burden, J. Praschifka, CD. Roberts, Phys. Rev. D 46 (1992) 2695. 
[14] P. Maris, Phys. Rev. D 52 (1995) 6087. 



17 



[15] M. Hermele, T. Senthil, M.P.A. Fisher, P.A. Lee, N. Nagaosa, X.-G. Wen, Phys. Rev. B 70 (2004) 
214437. 

[16] S. Takashima, I. Ichinose, T. Matsui, Phys. Rev. B 73 (2006) 075119. 

[17] Renormalization group analysis showed that the critical fermion flavor beyond which deconfinement 

occurs is as large as 36. See, F.S. Nogueira, H. Klcincrt, Phys. Rev. B 77 (2008) 045107. 
[18] Y. Nambu, G. Jona-Lasinio, Phys. Rev. 122 (1961) 345. 

[19] A.M. Polyakov, Gauge Fields and Strings (Harwood, Academic, Chur, Switzerland, 1987). 
[20] L.D. Landau, I.M. Khalatnikov, JETP 2 (1956) 69; B. Zumino, J. Math. Phys. 1 (1960) 1. 
[21] H. Gcorgi, E.H. Simmons, A.G. Cohen, Phys. Lett. B 236 (1990) 183; T. Kugo, M.G. Mitchard, Phys. 

Lett. B 282 (1992) 162; T. Ebihara, T. Iizuka, K.-I. Kondo, E. Tanaka, Nucl. Phys. B 434 (1995) 8; 

K.-I. Kondo, Phys. Rev. D 55 (1997) 7826. 
[22] G.-Z. Liu and G. Cheng, Phys. Rev. D 67 (2003) 065010; H. Jiang, G.-Z. Liu, G. Cheng, J. Phys. A 41 

(2008) 255402. 

[23] D.J. Lee, I.F. Hcrbut, Phys. Rev. B 66 (2002) 094512; O. Vafek, Z. Tesanovic, M. Franz, Phys. Rev. 

Lett. 89 (2002) 157003. 
[24] Y. Ran, M. Hermele, P.A. Lee, X.-G. Wen, Phys. Rev. Lett. 98 (2007) 117205. 
[25] W. Li, G.-Z. Liu, arXiv: 0907.2365vl. 
[26] C. Mudry, E. Fradkin, Phys. Rev. B 49 (1994) 5200. 
[27] G.-Z. Liu, Phys. Rev. B 71 (2005) 172501. 
[28] A. Durst, P. A. Lee, Phys. Rev. B 62 (2000) 1270. 
[29] V. P. Gusynin and V. A. Miransky, Eur. Phys. J. B 37 (2003) 363. 

[30] S. Yamashita, Y. Nakazawa, M. Oguni, Y. Oshima, H. Nojiri, Y. Shimizu, K. Miyagawa, K. Kanoda, 

Nature Phys. 4 (2008) 459. 
[31] M. Yamashita, N. Nakata, Y. Kasahara, T. Sasaki, N. Toneyama, N. Kobayashi, S. Fujimoto, T. 

Shibauchi, Y. Matsuda, Nature Phys. 5 (2008) 44. 



18 



